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Abstract 


A  single-sensor  hot-wire  probe  was  used  to  make  streamwise  measurements 
in  a  turbulent  plane  mixing  layer.  The  mixing  layer,  with  a  velocity  ratio  of 
approximately  2:1,  was  created  in  an  open  return  wind  tunnel  by  the 
insertion  of  a  curved  splitter  plate.  Conditions  were  studied  with  the  splitter 
plate  boundary  layers  tripped.  The  single-sensor  results  demonstrate  that 
this  mixing  layer  has  standard  velocity  field  statistical  properties.  The  data 
were  taken  at  several  Reynolds  numbers  in  the  fully  turbulent  flow  down¬ 
stream  of  the  mixing  point.  The  results  from  the  tripped  initial  bormdary 
layers /fully  turbulent  conditions  were  compared  with  the  temporally  de¬ 
veloping  mixing  layer  direct  numerical  simulation  results  of  Rogers  and 
Moser  (1994). 


Rogers,  M.  M.  and  R.  D.  Moser,  February  1994,  Direct  Simulation  of  a  Self-Similar  Turbulent  Mixing  Layer, 
Physics  of  Fluids  and  Fluid  Dynamics  6,  pp  903-923. 
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1.  Introduction 


The  purpose  of  this  research  was  to  obtain  a  thorough  knowledge  of  the 
structure  and  dynamics  of  the  streamwise  velocity  of  a  plane  mixing  layer 
using  a  single-sensor  hot-wire  probe,  developed  in  the  Turbulence  Re¬ 
search  Laboratory  (TRL)  at  the  University  of  Maryland  College  Park 
(UMCP).  The  results  of  the  measurements  were  compared  against  a  theo¬ 
retical  computation  of  a  similar  flow  by  Rogers  and  Moser  (1994). 

The  turbulent  plane  mixing  layer  is  technically  one  of  the  most  important 
types  of  generic  turbulent  flows.  It  is  the  basic  flow  field  created  to  en¬ 
hance  species  mixing  for  combustion  devices.  The  turbulent  mixing  layer 
is  also  the  imderlying  flow  in  many  types  of  fluid  instability  problems.  The 
results  of  this  project  could  lead  to  an  improvement  in  state-of-the-art  hir- 
bulence  models  used  in  computational  fluid  d)mamics  for  calculating  drag 
and  mixing. 


2.  Background 

Plane  mixing  layers  have  been  widely  studied,  both  experimentally 
(Bradshaw,  1966;  Wygnanski  and  Fiedler,  1970;  Miksad,  1972;  Batt,  1977) 
and  computationally  (Patnaik,  Sherman,  and  Corcos,  1976;  Riley  and 
Metcalf,  1980;  Jacobs  and  Pullin,  1989).  After  the  realization  by  Brown  and 
Roshko  (1974)  that  coherent  structures  dominate  some  of  the  mixing  dy¬ 
namics,  a  large  effort  was  imdertaken  to  investigate  these  structures  in  an 
effort  to  further  understand  the  onset  and  evolution  of  turbulence. 

Mixing  layers  play  an  important  part  in  many  engineering  applications: 
they  govern  the  rate  of  mixing  in  flow  reactors  and  combustion  chambers, 
and  they  are  the  dominant  cause  of  noise  from  modem  aircraft  propulsion 
systems.  Mixing  layers  are  also  of  interest  because  they  exhibit  self-similar 
behavior.  At  sufficiently  high  Reynolds  numbers,  the  mixing  layer  tends  to 
grow  linearly,  and  the  shapes  of  the  mean  velocity  and  turbulence  profiles 
are  independent  of  downstream  distance  when  scaled  with  appropriate  lo¬ 
cal  length  and  velocity  scales. 

The  region  between  the  parallel  streams  is  most  likely  the  simplest  free 
shear  flow  that  can  be  studied  since  the  driving  flows,  Uj^  and  LZj,  are  con¬ 
stant  throughout.  It  is  easily  reproduced  in  a  wind  tunnel  and  creates  large 
coherent  turbulent  stmctures  in  a  boundary-free  field.  A  large  body  of  ex¬ 
perimental  and  numerical  results  has  demoristrated  that  the  mixing  layer 
developed  by  the  meeting  of  two  streams  at  the  end  of  a  splitter  plate,  as 
depicted  in  figure  1,  is  dominated  by  mostly  two-dimensional  spanwise 
vortex  structures  rising  from  the  Kelvin-Helmholtz  instability  of  the  layer 
(Brown  and  Roshko,  1974;  Winant  and  Browand,  1974;  Corcos  and 
Sherman,  1984;  Metcalf,  Hussain,  Menon,  and  Hayakawa,  1987;  Huang 
and  Ho,  1990).  These  structures,  referred  to  as  "rollers,"  are  evenly  spaced 
by  some  length  I,  and  convect  downstream  at  a  constant  velocity,  but  grow 
in  size.  Eventually,  some  rollers  interact  by  rotating  about  each  other. 
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Figure  1.  Mixing  layer 
velocity  profile. 


These  interacting  rollers  will  pair,  forming  a  single  vortex  with  double  the 
original  spacing.  Pairing  occurs  repeatedly,  controlling  the  growth  rate  of 
the  mixing  layer. 


2.1  Experimental  Characteristics 

Brown  and  Roshko  (1974),  who  first  recognized  the  rollers,  recorded  their 
observations  with  a  series  of  shadowgraphs.  They  showed  that  the  roller 
persisted  down  the  length  of  their  test  region.  Since  the  mean  velocity  pro¬ 
files  seemed  self-preserving  and  the  spread  rate  was  linear.  Brown  and 
Roshko  suggested  that  this  was  the  unique,  fully  developed  form  of  a  tur¬ 
bulent  mixing  layer.  They  also  observed  that  some  rollers  would  merge,  or 
pair,  and  evolve  into  a  single  roller.  Winant  and  Browand  (1974)  discov¬ 
ered  that  at  low  Reynolds  numbers,  self-preservation  was  maintained  even 
after  several  stages  of  vortex  pairings,  without  complete  breakdown  to 
fully  turbulent  flow.  Dimotakis  and  Brown  (1976)  found  evidence  of  the 
roller  structures  at  high  Reynolds  numbers.  Chandrsuda,  Mehta,  Weir,  and 
Bradshaw  (1978)  speculated  that  the  rollers  in  the  mixing  layer  existed 
only  if  the  free-stream  turbulence  was  low.  They  rationalized  that  high  tur¬ 
bulence  in  the  free  stream  caused  the  rollers  to  break  down  early,  and  sig¬ 
nificantly  high  free-stream  turbulence  would  not  allow  the  rollers  to  form 
at  all.  They  also  suggested  that  the  rollers  were  not  a  universal  feature  of 
the  mixing  layer,  and  that  pairing  could  happen  along  part,  but  not  all,  of 
the  span  of  the  roller.  Browand  and  Troutt  (1980)  learned  that  the  rollers 
quickly  approached  an  asymptotic  state,  where  the  lateral  correlation 
length  became  proportional  to  the  local  mixing  length  (momentum  thick¬ 
ness).  They  also  observed  spanwise  irregularities,  but  argued  that  these  ir¬ 
regularities  were  a  product  of  pairing. 

Briedenthal  (1981)  visualized  the  development  of  the  rollers  with  an  aque¬ 
ous,  diffusion-limited  reaction.  He  observed  that  a  transition  occurred  in 
which  the  product  of  the  reactants  increased  by  an  order  of  magnitude.  He 
attributed  this  transition  to  the  development  of  small-scale,  three- 
dimensional  motions  in  the  flow.  Downstream  of  the  transition,  the  degree 
of  mixing  became  independent  of  the  local  Reynolds  number.  Heman  and 
Jiminez  (1982)  observed  that  most  of  the  entrainment  of  the  free-stream 
fluid  occurred  during  the  normal  life  of  the  rollers,  and  not  during  pairing. 
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Ho  and  Huang  (1982)  discovered  they  could  change  the  growth  rate  of  the 
rollers  with  subharmonic  perturbation  by  forcing  several  of  the  rollers  to 
merge  simultaneously.  This  agreed  with  Kelly  (1967),  who  proposed  that  a 
subharmonic  resonance  mechanism  could  feed  energy  to  large 
wavenumber  structures,  and  with  Pierrehumbert  and  Widnall  (1982),  who 
determined  that  an  array  of  rollers  were  unstable  to  subharmonic  distur¬ 
bances.  Jiminez  (1983)  concentrated  on  studying  some  streamwise  irregu¬ 
larities  that  appeared  on  the  rollers,  and  concluded  that  these  irregularities 
were  a  secondary  instability  and  actually  deformed  the  large  roller.  Lang 
(1985)  measured  the  instantaneous  values  of  the  streamwise  and  trans¬ 
verse  velocity  components,  spanwise  vorticity,  and  the  Re3molds  stress.  He 
discovered  that  the  spanwise  vorticity  component  had  values  opposite  to 
the  sense  of  the  mean  vorticity  most  of  the  time. 

Bernal  and  Roshko  (1986)  expanded  the  study  of  the  secondary  structures 
with  both  gas  and  water  facilities.  They  observed  the  same  type  of  streaks 
seen  earlier  by  Brown  and  Roshko  (1974).  By  visualizing  a  plane  normal  to 
the  streamwise  direction  with  the  use  of  laser-induced  fluorescence,  Bernal 
and  Roshko  visualized  the  secondary  instabilities  and  concluded  that 
these  instabihties  were  streamwise  vortices  superimposed  on  the  rollers.  A 
chemically  reacting  mixing  layer  was  used  by  Lasheras,  Cho,  and 
Maxworthy  (1986)  to  visualize  and  investigate  the  streamwise  vortices. 
They  determined  that  these  streamwise  vortices  were  a  product  of  up¬ 
stream  perturbations,  and  could  be  increased  in  intensity  by  further  per¬ 
turbing  the  upstream  flow.  The  streamwise  vortices  seemed  to  originate  in 
the  braid  region  between  the  rollers.  The  study  concluded  that  the 
streamwise  vortices  seemed  to  have  a  scale,  as  measured  by  the  size  of  the 
cores,  smaller  than  but  comparable  to  that  of  the  roller's  cores;  this  they 
significantly  contributed  to  the  entrainment  process  in  the  early  stages  of 
the  mixing  layer  development.  Lasheras  and  Choi  (1988)  performed  a 
similar  experiment,  except  they  introduced  a  sinusoidal  perturbation  by 
corrugating  the  splitter  plate  and  indenting  the  splitter  plate  edge.  They 
showed  that  the  mixing  layer  with  induced  streamwise  vorticity  grew  at  a 
greater  rate  than  the  flow  without  induced  streamwise  vorticity.  It  was 
again  surmised  that  the  streamwise  vortices  originated  within  the  braid  re¬ 
gion.  The  study  concluded  that  the  perturbed  vorticity  existing  within  the 
braids  underwent  axial  stretching  from  the  effect  of  the  strain  field  created 
by  the  evolving  rollers.  This  resulted  in  the  formation  of  vortex  tubes 
whose  axes  were  aligned  with  the  principle  direction  of  the  positive  strain 
field. 

Balint,  Wallace,  and  Vukoslavcevic  (1989)  demonstrated  the  ability  to 
measure  velocity  and  vorticity  directly  with  the  use  of  a  nine  hot-wire  ane¬ 
mometer  array.  They  were  able  to  collect  data  at  a  single  station  within  the 
shear  layer,  at  the  edge  of  the  layer,  and  in  the  free  stream.  The  dominant 
source  term  in  the  enstrophy  balance  was  due  to  vortex  stretching  from  the 
fluctuating  strain  rate. 

Bell  and  Mehta  (1990)  evaluated  the  development  of  the  mixing  layer 
originating  from  tripped  and  untripped  boundary  layers.  The  distance  to 
achieve  self-preservation  was  shorter  for  the  tripped  case,  while  the  mix- 
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ing  layer  growth  was  greater  for  the  untripped  case.  Bell  and  Mehta  attrib¬ 
uted  this  difference  to  the  presence  of  stream  wise  vortices  that  added  to 
the  entrainment  of  the  mixing  layer.  Bell,  Plesniak,  and  Mehta  (1992)  per¬ 
formed  a  spanwise  averaging  of  the  mixing  layer  and  found  large  varia¬ 
tions  in  the  mean  flow  and  turbulence  properties.  They  postulated  that  the 
variations  were  caused  by  the  presence  of  spatially  stationary  stream  wise 
vortices  that  ride  along  the  rollers.  They  cautioned  that  there  could  be  sig¬ 
nificant  misinterpretations  regarding  the  development  of  the  mixing  layer 
based  on  centerline  data  alone  and  proposed  that  subsequent  discrepan¬ 
cies  of  past  observances  regarding  growth  rates  and  turbulence  properties 
could  be  explained  by  the  lack  of  spanwise  averaging. 

Numerical  Simulations 

Corcos  and  Sherman  (1984)  simulated  the  two-dimensional  evolution  of 
the  roller  using  pseudo-spectral  methods,  and  were  able  to  simulate  pair¬ 
ings.  Corcos  and  Lin  (1984)  showed  that  the  rollers  were  unstable  due  to 
three-dimensional  disturbances  that  resulted  in  their  deformation  and  the 
formation  of  the  streamwise  vortices  in  the  braid  region.  Ashurst  and 
Meiburg  (1988)  simulated  the  mixing  layer  using  a  discrete  vortex  method. 
They  included  two  signs  of  vorticity  to  account  for  the  weak  boundary 
layer  leaving  the  splitter  plate.  The  result  was  a  simulation  in  closer  agree¬ 
ment  with  experimental  observations.  Ashurst  and  Meiburg  were  able  to 
show  the  generation  of  streamwise  vortices  originating  within  the  braid 
region.  Ghoniem  and  Heidarinejad  (1991)  found  that  there  was  a  strong 
similarity  between  the  distribution  of  product  concentrate  and  vorticity  at 
a  wide  range  of  Damkohler  numbers  (the  Damkohler  number  is  the  ratio 
between  convective  flow  and  chemical  reaction  times).  They  attributed  the 
similarity  to  the  entrainment  associated  with  pairing. 

Moser  and  Rogers  (1991, 1992, 1993)  performed  several  direct  numerical 
simulations  with  a  variety  of  initial  conditions.  They  observed  that  pairing 
the  rollers  caused  them  to  undergo  a  transition  to  small-scale  turbulence. 
Pairing  inhibited  the  growth  of  infinitesimally  small  disturbances  and  trig¬ 
gered  the  transition  to  turbulence  in  highly  three-dimensional  flow.  They 
concluded  that  streamwise  vortices  were  part  of  the  initial  cascade  to 
small-scale  turbulence.  They  also  found  that,  in  the  self-preserving  regime, 
the  streamwise  vortices  were  absent,  and  the  mixing  layer  growth  was 
reduced. 

Rogers  and  Moser  (1992)  studied  in  great  detail  the  origination  and  evolu¬ 
tion  of  the  streamwise  vortices,  also  referred  to  as  "rib"  vortices.  They 
foimd  that  the  spanwise  structures  evolved  into  a  corrugated  roller,  which 
was  consistent  with  the  observations  of  Jiminez  (1983),  with  vortex  stretch¬ 
ing  creating  strong  spanwise  vorticity  in  a  cup-shaped  region  at  the  bends 
of  the  rollers.  For  strong  initial  perturbations,  the  rib  vortices  collapsed 
into  tight  axisymmetric  vortices.  The  rib  vortex  lines  connected  to  other 
nearby  rib  vortices,  and  were  kinked  in  the  direction  opposite  to  the  roller 
vortex  direction.  Confirming  earlier  observations  of  Metcalf,  Orszag, 
Brachet,  Menon,  and  Riley  (1987),  and  Huang  and  Ho  (1990),  Rogers  and 


Moser  found  that  pairing  prevented  the  growth  of  three-dimensionality  by 
halting  oversaturation. 

Using  a  turbulent  boundary  layer  simulation  (Spalart,  1988)  as  the  mixing 
layer  initial  condition,  Rogers  and  Moser  (1994)  performed  three  simula¬ 
tions,  one  without  any  disturbance  in  the  initial  conditions,  and  the  other 
two  with  varying  disturbances  in  the  initial  conditions.  They  found  that, 
after  a  development  stage,  the  case  with  weak  and  no  initial  disturbance 
evolved  in  a  self-preserving  manner.  This  self-similar  period  was  charac¬ 
terized  by  the  absence  of  roller  pairings,  a  lack  of  rib  vortices  in  the  braid 
region,  and  scalar  mixing  characterized  by  "marching"  probability  density 
functions  (PDFs)  (the  largest  probable  value  varies  across  the  mixing  layer, 
with  the  largest  probable  value  on  each  side  of  the  mixing  layer  being 
closer  to  the  free-stream  value  of  that  side).  The  simulation  using  the 
strong  initial  disturbance  approached  self-preservation,  but  demonstrated 
strong  organized  pairings  of  the  rollers,  clearly  defined  regions  within  the 
braids  of  strong  streamwise  vorticity,  and  scalar  PDFs  that  were 
"nonmarching."  The  statistics  and  structures  from  several  experiments 
were  in  better  agreement  with  the  results  of  the  simulations  that  did  not 
exhibit  organized  pairings. 


3.  Experiment  Program 

Clearly,  research  on  planar  mixing  layers  has  been  quite  extensive  for  over 
20  years.  Fortunately,  state-of-the-art  development  in  hot-wire  anemom- 
etry  has  more  recently  allowed  reasonably  well-resolved  instantaneous 
measurements  of  the  velocity  vector,  vorticity  vector,  and  strain  rate  tensor 
(Vukoslavcevic,  Wallace,  and  Balint,  1991).  TTiis  development  provided  the 
opportunity  to  obtain  new  information  regarding  aspects  of  turbulent  phe¬ 
nomena  within  the  mixing  layer.  To  date,  much  of  the  experimentally 
measured  data  have  concentrated  on  the  streamwise  and  transverse  veloci¬ 
ties,  the  associated  stresses,  and  the  spanwise  vorticity.  Except  for  visual 
observations,  there  was  little  qualitative  information  and  very  few  quanti¬ 
tative  measurements  of  the  secondary  vortical  structures,  much  less  any 
statistical  information  related  to  their  origin  and  evolution.  There  is  a  data¬ 
base  of  direct  numerical  simulations  that  provides  considerable  informa¬ 
tion  (see  work  by  Rogers  and  Moser).  The  conclusions  drawn  from  these 
numerical  simulations  needed  to  be  verified  experimentally  and,  if  pos¬ 
sible,  extended. 

I  modified  the  wind  tunnel  m  the  TRL  at  UMCP  with  a  special  insert,  simi¬ 
lar  in  design  to  the  one  used  successfully  by  Dimotakis  et  al.  (1976)  and 
Lang  (1985).  The  insert  separated  the  conditioned  flow  coming  from  the 
wind  tunnel  contraction,  accelerating  the  flow  on  one  side  while  decelerat¬ 
ing  it  on  the  other.  The  end  of  the  insert  provided  a  smooth  meeting  of  the 
two  laminar  flows,  with  a  velocity  ratio  of  approximately  2:1,  and  created  a 
mixing  layer. 
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A  series  of  single-sensor  anemometer  measurements  and  flow  visualiza¬ 
tions  were  made  to  determine  if  the  turbulence  characteristics  of  the  mix¬ 
ing  layer  were  consistent  with  the  literature.  Specifically  the  goal  was  to 
match  as  closely  as  possible  the  simulated  conditions  of  Rogers  and  Moser 
(1994).  To  evaluate  the  characteristics  of  the  wind  tunnel  facility  a  series  of 
measurements  were  made  with  a  single  hot-wire  anemometer  traversing 
the  mixing  layer.  The  information  obtained  from  the  sensor  was  resolved 
to  determine  the  turbulence  statistical  characteristics,  including  moments 
and  one-dimensional  spectra. 

3.1  Preliminary  Characterization 

Following  the  method  described  by  Dimotakis  et  al.  (1976)  and  Lang 
(1985),  and  used  in  the  GALCIT  Facility,  an  insert  for  the  wind  tunnel  at 
the  UMCP's  TRL  was  designed,  fabricated,  and  installed.  I  obtained  a  se¬ 
ries  of  single  hot-wire  measurements  to  evaluate  the  ability  of  the  insert/ 
wind-tunnel  facility  to  generate  a  planar  mixing  layer. 

3.1.1  Initial  Computational  Analysis 

To  optimize  the  design  of  the  insert,  I  performed  a  numerical  computation 
of  a  two-dimensional,  steady  flow  using  a  commercial  software  called 
ALGOR™  (1992).  I  believed  this  would  give  a  reasonable  approximation  of 
the  flow  to  be  encountered.  The  model  consisted  of  3746  elements  and  was 
solved  on  a  Pentium-processor  PC  in  about  12  hours  using  the  Penalty 
Parameter  Method  (ALGOR,  1992).  I  used  a  uniform  flow  of  1  m/s  as  the 
input  boundary  condition  on  the  left  side,  and  enforced  the  no-slip  condi¬ 
tion  along  the  top  and  bottom  surfaces  of  the  insert.  I  used  a  rounded  lead¬ 
ing  edge  for  the  calculation,  which  helped  fix  the  location  of  the  leading 
edge  stagnation  point  and  reduced  the  chance  of  separation.  The  calcula¬ 
tion  was  a  pseudo-steady-state  calculation;  therefore  the  rollers  would  not 
appear  in  the  results.  Figure  2  shows  the  model's  mesh,  the  results  of  the 
calculations  of  velocity  magnitude,  streamwise  and  transverse  velocity 
components,  vorticity,  and  pressure.  The  flow  is  represented  moving  from 
left  to  right.  The  color  bar  on  the  right  of  each  plot  shows  the  spatial  varia¬ 
tion  of  each  flow  component. 

Figure  2f  shows  that  an  unequal  pressure  distribution  exists  along  the  en¬ 
trance  plane.  The  pressure  drops  a  great  deal  as  it  passes  below  the  insert. 
Above  the  insert,  the  pressure  does  not  drop  as  much,  similar  to  the  pres¬ 
sure  distribution  at  a  contraction.  This  difference  in  pressure  distribution 
drives  the  incoming  fluid  down  into  the  lower  pressure  region.  This  pro¬ 
duces  a  strong  negative  vertical  component  to  the  velocity  vector,  as  seen 
in  figure  2d.  The  vorticity  in  figure  2e  also  indicates  that  an  area  of  recircu¬ 
lation  exists  at  the  lower  insert  entrance.  The  change  in  the  streamwise  ve¬ 
locity  component  in  the  transverse  direction  is  very  small,  away  from  the 
insert  and  walls.  A  transverse  change  in  the  streamwise  velocity  compo- 
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Figure  2.  Two- 
dimensional 
calculations  of  flow 
around  insert  in 
tunnel:  (a)  model 
mesh,  (b)  velocity 
magnitude, 

(c)  streamwise 
velocity  component, 

(d)  transverse  velocity 
component, 

(e)  spanwise  vorticity, 
and  (f)  pressure. 


nent  results  in  changes  in  the  spanwise  vorticity.  The  largest  changes  in 
vorticity  occur  at  the  leading  edge  of  the  insert  and  along  the  convex 
curved  section.  The  resulting  angle  of  attack  of  the  flow  around  the  leading 
edge  forms  a  separation  bubble  at  the  low  side  of  the  leading  edge.  As  a 
result,  the  insert  acts  as  a  bluff  body,  with  a  weak  separation  bubble  filling 
the  concave  region  of  the  expansion.  The  streamwise  velocity  is  acceler¬ 
ated,  as  shown  in  figure  2c,  as  the  fluid  passes  the  insert.  There  is  little  dif¬ 
ference  in  the  free-stream  velocities  downstream  of  the  splitter  plate. 

The  simulation  exposes  a  serious  deficiency  in  the  insert  design.  Fortu¬ 
nately,  Dimotakis  et  al.  (1976)  and  Lang  (1985),  who  encountered  the  same 
problem,  eliminated  the  separation  at  the  leading  edge  by  placing  a  flow 
restriction  in  the  expansion  side.  In  this  manner,  they  were  able  to  match 
inlet  and  outlet  pressures.  They  placed  perforated  plates  with  50-percent 
blockage  on  the  expansion  side,  where  separation  might  otherwise  occur. 
The  plates  were  curved  to  approximately  follow  the  shape  of  a  potential 
flow  line  of  an  inviscid  flow  through  the  expansion.  Fine  adjustments  of 
the  pressure  drop  were  accomplished  with  an  identical  perforated  plate 
placed  over  the  first.  Sliding  one  plate  vertically  relative  to  the  other  varied 
the  percent  blockage.  This  caused  a  small  problem  in  that  the  flow  direc¬ 
tion  was  slightly  changed,  and  this  had  to  be  accounted  for  by  tilting  the 
plate.  Lang  (1985)  foimd  that  by  reducing  the  blockage  near  the  insert,  the 
boundary  layer  could  be  made  smaller  at  the  end  of  the  expansion. 


3.1.2  Final  Design  of  Insert 

An  insert  was  fabricated  from  two  0.00318-m  (1/8-in.)  thick  sheets  of 
plexiglass,  each  1.22  x  0.91  m.  A  semi-circular  leading  edge  and  tapering 
trailing  edge  were  used  to  connect  the  sheets,  as  shown  in  figure  3a.  Six 
0.00318-m  diam.  steel  rods  passed  through  the  gap  in  the  plexiglass  sheets 
and  secured  the  insert  to  the  wind  tunnel  walls.  The  insert  was  placed  in 
the  upstream  part  of  the  wind  tunnel,  about  2  m  downstream  of  the  end  of 
the  contraction.  As  shown  in  figure  3b,  the  insert  served  to  separate  the 


Figure  3.  UMCP 
wind  tunnel  insert; 
(a)  components 
and  (b)  placement. 


(b) 


8 


flow,  accelerating  the  flow  on  one  side  and  decelerating  it  on  the  other.  The 
downstream  part  of  the  insert  served  as  the  splitter  plate,  and  tapered  to  a 
sharp  point  with  an  angle  of  3.80°  on  either  side.  The  insert  splits  the  flow 
such  that  there  is  a  velocity  ratio  of  about  2:1  at  the  splitter  plate  exit  plane. 

It  was  clear  that  a  flow  restriction  on  the  deceleration  side  was  necessary  to 
equalize  the  pressure  at  the  insert  entrance  plane.  As  shown  in  figure  3b,  a 
screen  of  about  38  percent  blockage  was  placed  at  the  most  likely  point  of 
separation  on  the  deceleration  side  for  two  reasons:  to  equalize  the  pres¬ 
sure  differential  across  the  insert  and  to  prevent  flow  separation  in  the  re¬ 
gion  of  adverse  pressure  gradient. 

3.1.3  Flow  Visualization 

Smoke  was  used  to  visualize  the  flow,  to  help  determine  if  the  insert  was 
effective.  We  created  smoke  by  burning  a  smoke-generating,  non-toxic 
incense  material  in  a  combustion  chamber.  We  pressurized  the  chamber  to 
70  kPa  and  extended  a  hose  from  the  top  of  the  combustion  chamber  to  a 
manifold  via  a  heat  exchanger.  The  manifold  distributed  the  smoke  to  four 
smaller  lines,  which  emptied  into  a  plenum  inside  the  wind  tunnel  insert. 
We  formed  the  insert  plenum  from  the  gap  left  by  the  plexiglass  sheet,  as 
described  in  section  3.1.2.  The  0.00318-m  gap  was  approximately  0.4  m 
wide  and  1.22  m  deep,  creating  a  plenum  volume  of  1.55  liters.  About 
0.4  m  upstream  of  the  trailing  edge,  cut  into  the  high  speed  side  of  the  in¬ 
sert,  is  a  0.00159-m-wide,  0.2-m-long  slit,  where  the  smoke  seeps  into  the 
flow  at  about  0.1  m/s.  A  flash  bulb,  placed  downstream  but  facing  up¬ 
stream,  illuminated  the  smoke  as  it  was  convected  downstream.  Figure  4 
shows  a  series  of  such  captured  images.  In  each  case,  the  wind  tunnel  was 
running  at  different  speeds,  indicated  by  the  blowers  speed  in  revolutions 
per  minute  (rpm). 

The  images  were  captured  on  200- ASA  film,  with  the  /-stop  set  to  2.5  and 
the  shutter  speed  set  to  30.  This  is  a  relatively  slow  shutter  speed,  but  with 
the  low  ASA  film,  this  was  the  fastest  shutter  speed  it  was  possible  to  use 
without  losing  the  image.  As  a  result,  the  higher  velocity  flows  (fig.  4d 
and  e)  are  somewhat  smeared.  Despite  these  limitations,  spanwise  vortical 
concentrations  are  still  visible  at  the  wind  tunnel's  highest  speed. 

3.2  Hot-Wire  Measurements 

Hot-wire  anemometry  is  still  the  method  used  most  often  to  measure  the 
fluctuating  velocity  components  of  fluids.  The  high  spatial  resolution,  fast 
response  time,  and  low  fabrication  cost  make  hot  wire  sensors  ideal  for  this 
investigation  (Blackwelder,  1981).  A  series  of  single  hot-wire  measure¬ 
ments  were  made  to  characterize  the  flow  in  the  wind  tunnel,  and  then 
compared  to  the  velocity  characteristics  of  the  results  of  Rogers  and  Moser 
(1994).  It  was  discovered  that  the  contraction  side  of  the  splitter  plate  re¬ 
quired  a  trip  to  produce  a  turbulent  boundary  layer.  Once  established,  the 
resulting  velocity  measurements  were  very  close  to  the  numerical  results, 
as  well  as  the  results  of  Bell  and  Mehta  (1990). 
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Figure  4.  Streaklines 
from  smoke  shown 
using  flashbulb  with 
200-ASA  film:  (a)  200 
rpm,  (b)  270  rpm, 

(c)  400  rpm,  (d)  500 
rpm,  and  (e)  700  rpm. 


3.2.1  Measurement  Results 

A  single  hot  wire  was  used  to  obtain  instantaneous  velocity  measurements 
for  the  condition  of  upstream  velocity  =  2.46  m/ s.  The  measurements 
were  made  at  0.005-m  increments  along  a  transverse  line  ranging  from 
0.15  m  below  (slow  side)  to  0.15  m  above  (fast  side)  the  trailing  edge.  A 
complete  set  of  measurements  was  obtained  from  just  at  the  downstream 
trailing  edge  of  the  splitter  plate  (Xj  ~  0  m)  to  1.1  m  downstream,  at  0.1-m 
intervals.  The  hot  wire  was  moved  to  each  vertical  location  with  a  traverse 
mechanism,  and  held  for  60  s  while  the  acquisition  system  sampled  at 
1  kHz.  The  data  were  analyzed  for  mean  velocity,  variance,  standard  de¬ 
viation,  skewness,  and  kurtosis.  Figures  5a  through  e  are  the  resulting 
plots. 
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(a)  mean  velocity, 

(b)  dUj^/dX2r 

(c)  variance. 


Figure  5.  Contour 
results  of  single  hot¬ 
wire  measurements: 

(d)  skewness,  and 

(e)  kurtosis. 


These  figures  consist  of  a  three-dimensional  data  point  plot,  followed  by  a 
colored  contour  plot.  The  data  points  are  plotted  in  a  perspective  view 
with  the  measured  quantity  represented  as  vertical  height.  The  color  bar 
on  the  right  of  each  contour  plot  shows  the  spatial  variation  of  each  flow 
component.  They  show  that  the  mixing  layer  grows  linearly  from  about 
0.3  m  downstream  of  the  splitter  plate.  Table  1  lists  some  resulting  param¬ 
eters.  These  results  are  discussed  later. 

There  are  two  convenient  lengths  for  scaling  that  can  be  derived  from  the 
mean  velocity  distribution.  The  vorticity  thickness  5^  is  defined  as 


dU, 


9xi 


max 


(1) 


It  is  an  appropriate  measure  of  the  growth  of  the  turbulent  mixing  layer. 
The  momentum  thickness  0is  defined  in  the  usual  way,  as 


J—oo 


UizHc] 
AU 


\dXn 
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Table  1.  Characterization  results. 


Description 

Variable 

Numerical  value 

Free-stream  velocity,  fast 

3.19  m/s 

Free-stream  velocity,  slow 

1.69  m/s 

Velocity  ratio 

Ui/U, 

0.53 

Convection  velocity 

2.44  m/s 

Momentum  thickness  (m) 

d 

-0.000504  m  +  0.0165  x  1 

Reynolds  number  (based  on  6) 

Re^ 

80.7  + 1655.1 /mxl 

Vorticity  thickness  (m) 

Sco 

0.005484  m  +  0.06545  x  1 

Reynolds  number  (based  on  5^ 

427.99  +  6523.5/m  x  1 

Nondimensional  growth  rate  of 

0.0531 

Strouhal  number 

St 

0.12 

Wavelength 

A 

0.133  m 

Maximum  turbulence  level 

/AD  X 100% 

max 

27.8  % 

where  J  0.5(U;j  +  U^),  and  represents  the  flow  of  momentum  through  the 
mixing  layer  compared  to  that  which  would  occur  if  no  shear  layer  was 
present.  Since  the  measurements  were  taken  at  0.005-m  intervals  in  the  X2 
direction,  the  momentum  thickness  was  calculated  discretely  as 


0=1 

i 


1 

4 


AU  j 


2' 


0.005  m  , 


where  i  is  the  number  of  points  and  associated  location  along  X2.  Each 
length  scale  can  be  used  in  a  Re3molds  number.  Two  such  Re)molds  num¬ 
bers  are  defined:  the  Reynolds  number  based  on  vorticity  thickness. 


Re 


(4) 


and  the  Reynolds  number  based  on  momentum  thickness. 


(5) 

The  momentum  thickness  is  often  used  for  self-similar  scaling  rather  than 
vorticity  thickness  because  momentum  thickness  is  less  sensitive  to  statis¬ 
tical  noise.  Other  investigators  have  used  the  vorticity  thickness  to 
nondimensionally  scale  the  vertical  distance  across  the  shear  layer.  The  ra¬ 
tionale  is  that  this  scaling  takes  into  account  some  of  the  deviation  from 
linear  growth  of  the  shear  layer  near  the  trailing  edge.  Figure  6  shows 
growth  of  the  vorticity  thickness  with  distance  downstream,  and  the  corre¬ 
sponding  growth  of  the  Reynolds  numbers  based  on  5^.  Each  plot  also  has 
a  least-squares  linear  fit  to  the  data. 

The  growth  of  the  vorticity  thickness  and  the  related  Reynolds  number  are 
linear.  The  least-squares  linear  fit  of  the  data  yields  a  growth  rate,  d6^/dx2 
(the  slope),  and  the  initial  value  (intercept).  The  least-squares  linear  fit  to 
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Figure  6.  Growth  of 

(a)  vorticity  thickness 
with  thickness 
downstream  and 

(b)  Reynolds  number. 


0  0.2  0.4  0.6  0.8  1  12 

Downstream  location  (m) 


Downstream  location  (m) 


1.2 


the  data  set  appear  in  table  1.  The  nondimensional  growth  rate  of  the  vor¬ 
ticity  thickness  (r^  is  defined  as 


^^~AUdX2  ■ 


(6) 


Using  the  least-squares  fit  for  the  slope  (0.06545  m/m),  r^-  0.1065.  This  is 
a  little  higher  than  the  value  of  0.098  from  Wygnanski  and  Fiedler  (1970). 
Brown  and  Roshko  (1974)  attained  a  value  of  0.0905.  Their  lower  value 
may  be  due  to  the  fact  that  Brown  and  Roshko's  Reynolds  number,  based 
on  vorticity  thickness  (=342000  Xj),  was  much  higher  due  to  their  increased 
convective  velocity  (factor  of  2)  and  density  (factor  of  7).  The  simulation  by 
Rogers  and  Moser  (1994)  resulted  in  =  0.0672.  Similar  plots  were  made 
for  the  momentum  thickness,  6,  along  with  its  related  Reynolds  number 
(see  fig.  7). 


14 


Figure  7.  Growth  of 

(a)  momentum 
thickness  with 
thickness 
downstream  and 

(b)  associated 
Re5molds  number. 
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From  the  trailing  edge  out  to  about  0.3  m,  the  wake  defect  caused  by  the 
screen  on  the  expansion  produces  mean  velocity  profiles  that  are  not  simi¬ 
lar  to  the  profiles  further  downstream,  as  can  be  seen  in  figure  5a.  In  fact, 
the  profiles  behave  more  like  the  combination  of  two  flat-plate  boundary 
layers  and  the  wake  from  a  bluff  body.  The  momentum  thickness  may  be 
more  appropriately  calculated  (and  is  here)  out  to  0.2  m  downstream  of  the 
splitter  plate.  Self-similarity  is  attained  by  0.3  m  downstream,  and  the  mo¬ 
mentum  thickness  grows  as  expected.  Kim  (1989)  shows  a  linearly  grow¬ 
ing  momentum  thickness  starting  at  0.3  m.  Winant  and  Browand  (1974) 
show  the  momentum  thickness  evolving  from  the  trailing  edge.  This  is  eas¬ 
ily  attributed  to  the  different  method  employed  in  creating  the  mixing 
layer  from  two  separate  flows,  and  exploiting  the  lower  kinematic  viscos¬ 
ity  of  water.  Wmant  et  al.  explained  that  0  initially  grows  as  the  square  root 
of  the  downstream  distance.  Where  the  mixing  layer  transitions  to  turbu¬ 
lence,  6  grows  linearly.  The  conditions  of  Rogers  and  Moser  (1994)  was  an 
initial  Re^  and  Reg  of  1370  and  800,  growing  to  10,800  and  2,420, 
respectively. 
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Figure  8.  Normalized 
mean  velocity 
profiles. 


Figure  9.  Normalized 
Reynolds  normal 
stress  in  streamwise 
direction. 


An  important  test  for  performance  of  the  facility  to  generate  a  mixing  layer 
matching  the  numerical  conditions  of  Rogers  and  Moser  (1994)  is  whether 
the  data  collapses  in  a  self-preserving  manner.  Rogers  and  Moser  collapsed 
the  mean  velocity  data  by  plotting  the  nondimensional  mean  streamwise 
velocity  (U^  -  U^)/AU,  plotted  against  the  nondimensional  vertical  dis¬ 
tance,  \  =  x^lQ,  where  AU=U^-Ui.  As  shown  in  figure  8,  the  mean  veloc¬ 
ity  from  the  preliminary  study  shows  self-similarity  from  about  0.3  m 
downstream  of  the  splitter  plate.  The  collapse  of  the  data  is  in  reasonable 
agreement  with  the  results  of  Rogers  and  Moser. 

Rogers  and  Moser  indicate  that  a  more  sensitive  indicator  of  self-similarity 
is  the  collapse  of  the  Reynolds  stress  profiles  ww  ^  with  The  appropriate 
scale  used  was  AU^.  Their  data  showed  a  well-collapsed  distribution  of  the 
streamwise  Reynolds  normal  stress  at  five  different  locations  within  their 
self-similar  regime.  Figure  9  shows  the  normalized  data  results  from  this 
preliminary  study.  Much  like  the  results  of  Roger  and  Moser,  the  data  at  0.9 
through  1.1  m  very  nearly  collapse,  closely  matchiag  the  results  of  Rogers 
and  Moser's  peak  self-similar  value  of  wu  j  / AU^,  going  to  0  at  x  =  ±4. 
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3.2.2  Initial  Conditions 


An  initial  profile  of  the  mean  streamwise  flow  0.002  m  downstream  of  the 
splitter  plate  is  shown  in  figure  10.  One  can  see  that  the  mean  flow  exhibits 
a  larger  boundary  layer  profile  on  the  expansion  side,  and  the  variance 
seems  to  grow  with  the  velocity  defect.  It  is  theorized  that  there  is  some 
separation  along  the  concave  curvature  of  the  splitter  plate.  On  the  con¬ 
traction  side,  the  boundary  layer  is  initially  laminar  and  small. 

The  streamwise  mean  velocity  and  variance  data  of  the  expansion  side  of 
the  splitter  plate  were  plotted  to  evaluate  the  similarity  to  the  expected 
logarithmic  form  given  by  Clauser  (1956)  fit  to  a  turbulent  boundary  layer 
on  a  smooth,  solid  surface.  The  results  are  shown  in  figure  11,  where  the 
velocity  is  normalized  by  the  friction  velocity  u^,  defined  as 


where  the  velocity  and  length  are  normalized  as 


X2Ur 

Comparison  of  the  resulting  fit  indicates  the  mean  streamwise  velocity  and 
the  variance  with  Balint,  Wallace,  and  Vukoslavcevic  (1991)  are  not  behav¬ 
ing  exactly  as  a  turbulent  boundary  layer  on  a  flat  plate.  There  is  obviously 
a  contribution  from  the  relaxing  pressure  gradient  from  the  splitter  plates 
curvature.  It  is  expected  that  the  velocities  would  not  collapse  as  neatly  as 
the  flat  plate  family,  because  of  an  existing,  relaxing  pressure  gradient. 

Clauser  (1956)  named  a  similar  class  of  turbulent  boundary  layer  flows 
with  pressure  gradients  as  equilibrium  pressure  gradient  flows.  These  were 


Figure  10.  Initial 
streamwise  profiles. 


X2  (m) 
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Figure  11.  Boundary 
layer  profile  of 
expansion  surface: 

(a)  mean  streamwise 
velocity  and 

(b)  variance  profile. 


used  to  generalize  turbulent  boundary  layer  flows  around  curved  surfaces, 
much  like  the  splitter  plate.  Coles  (1956)  extended  the  law  of  the  wake, 

M  +  =  5.61ogy  +  +  4.9  +  5w(^j,  (9) 

to  equilibrium  pressure  gradient  by  letting  the  wake  parameter,  P,  be  a  func¬ 
tion  of  the  displacement  thickness,  wall  shear,  and  the  pressure  gradient: 


n  =  / 


dP] 

dX2j  ■ 


(10) 


Here,  <5  is  the  boundary  layer  thickness;  is  the  displacement  thickness  and 
is  defined  as 


while  the  wall  shear  is  defined  as 


9^2 


(11a) 


(11b) 


White  (1974)  proposed  that 


n  =  0.8 


^  dP 
dx2 
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(12) 


Letting  k=  0.41  (Schetz,  1984),  Cole  (1956)  proposed  that 

resulting  in  the  law  of  the  wake  as 


(13) 


u  ■'■  =  5.6  log  y  ■^  +  4.9  +  3.9024 


'iL  dP] 
dX2j 


(14) 


This  law  states  that  a  pressure  gradient  in  the  X2  direction  greater  in  magni¬ 
tude  than  that  of  a  flat  plat  will  result  in  larger  values  of  «■*■.  The  magnitude 
of  the  pressure  gradient  determines  the  degree  of  deviation  from  the 
smooth  flat  plate  turbulent  boundary  layer  profile.  In  this  instance,  on  both 
sides  of  the  splitter  plate,  the  pressure  gradient  is  relaxing  with  increasing 
Xj.  As  a  result,  the  added  value  of  the  pressure  gradient  contributes  to  the 
small  deviation  from  the  Clauser  fit  (Clauser,  1956),  as  one  can  observe 
from  figure  11b. 

The  variance  of  the  vmtripped  flow  is  already  in  excess  of  the  expected  val¬ 
ues  from  Balint  et  al.  (1991).  Tripping  the  expansion  side  would  create  even 
greater  variance  and  turbulent  kinetic  energy.  A  typical  peak  variance  of  a 
turbulent  boimdary  layer  is  about  u'/u^=  3.12  at  ~  11.4.  In  this  case,  is 
estimated  to  be  about  0.097  m/s. 


This  degree  of  turbulence  originating  from  the  expansion  side  of  the  split¬ 
ter  plate  eliminated  the  possibility  of  performing  ^e  simulation  with  lami¬ 
nar  initial  conditions.  The  contraction  side  was  tripped  by  placing  a 
0.01-m-diameter  rod  on  the  contraction  surface  about  0.1  m  downstream  of 
the  splitter  plate  leading  edge.  The  single  sensor  was  used  to  measure  the 
streamwise  velocity  on  the  contraction  side  about  0.005  m  upstream  of  the 
trailing  edge.  Figure  12  shows  the  measurements  obtained. 


The  boimdary  layer  profile  of  the  contraction  side  compares  much  more 
favorably  with  the  results  of  Balint  et  al.  (1991).  This  is  because  the  relax¬ 
ation  of  the  pressure  gradient  is  not  as  pronounced.  On  this  side,  a  typical 
peak  variance  of  a  turbulent  boundary  layer  is  about  «'/«(  =  2.37  at  y*"  = 
12.4,  and  u^is  estimated  to  be  about  0.19  m/s.  Table  2  lists  the  boimdary 
layer  properties  for  each  side.  Each  property  is  defined  as 


dU, 


dXn 


X2  =  0 


V 


C.  =  2- 


U 


(15) 


where  Cy  is  the  coefficient  of  drag. 
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Figure  12.  Boundary 
layer  profile  of 
tripped  contraction: 

(a)  mean  streamwise 
velocity  and 

(b)  variance  profile. 


Table  2.  Turbulent 
boundary  layer. 

Flow  side 

“t 

(m/s) 

dUi 

^^2  ^2  =  0 

C 

F 

u[ 

“r 

max 

«T 

Expansion 

Contraction 

0.097 

0.19 

611.37  s-’ 
2344.7  s-’ 

0.005958 
0.005706 ' 

3.12 

2.37 

0.0546 

0.0534 

The  effect  of  the  pressure  gradient  is  small.  As  a  result,  I  decided  to  per¬ 
form  simulations  with  the  contraction  side  tripped  and  the  expansion  side 
untripped.  This  allowed  me  to  use  our  existing  facilities  to  closely  emulate 
the  initial  condition  of  a  turbulent  boundary  layer  matching  the  simulation 
of  Rogers  and  Moser  (1994). 

3.3  Roller  Passage 

Most  of  the  references  discuss  the  frequency  of  the  passage  of  the  rollers. 
Kim  (1989)  showed  the  spectral  content  of  the  streamwise  and  transverse 
velocity  components.  His  plots  exhibited  a  smooth  cresting  between  30 
and  40  Hz  for  the  streamwise  velocity  spectrum,  while  there  was  a  defini¬ 
tive  peak  at  about  40  Hz  for  the  transverse  velocity  component. 

As  seen  in  figure  4,  there  is  strong  visual  evidence  that  rollers  are  produced 
in  this  facility.  Figure  13a  shows  a  sample  of  the  fluctuating  data,  taken  at 
the  edge  of  the  faster  part  of  the  shear  layer,  at  =  0.2  m. 
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Figure  13.  Sample 
data  and  normalized 
spectrum:  (a)  sample 
of  fluctuating  velocity 
data  from  single  hot¬ 
wire  anemometer  and 
(b)  normalized  power 
spectrum  of  velocity 
data. 

0.0  0.2  0.4  0.6  0.8  1.0 

Time  (s) 
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There  is  a  dominant  frequency  from  the  coherent  spanwise  vorticity  con¬ 
tinued  in  the  rollers.  If  the  frequency  of  passage  of  the  rollers  is  directly  re¬ 
lated  to  the  downstream  distance,  one  could  collapse  the  spectral  plots  by 
use  of  the  mixing  layer  scale,  6  or  5^.  Figure  13b  shows  the  spectral  energy 
content  of  this  data.  Note  the  apparent  collapse  of  the  spectral  profiles 
taken  from  from  0.3  to  1.1  m.  This  is  similar  to  the  result  of  Hussain  and 
Zaman  (1985).  Typically,  the  dominant  periodicity  of  the  passage  of  the 
rollers, /q,  is  expressed  nondimensionally  as  the  Strouhal  number 


St  = 


(15) 


The  resulting  Strouhal  number,  using  9  as  the  length  scale,  is  St  =  0.12  at  0.3 
to  1.1  m  downstream  of  the  trailing  edge.  Since  the  rollers  convect  down¬ 
stream  at  an  approximately  constant  velocity,  St  will  remain  constant.  The 
St  does  vary  from  x-^  =  0.3  m  to  1.1  m,  suggesting  evidence  of  pairing;  the 
rollers  grow  linearly,  yet  convect  downstream  at  a  constant  velocity,  even¬ 
tually  crowding,  and  then  pairing  occurs.  This  mechanism  is  what  allows  A 
to  remain  constant  in  the  self-similar  region  of  the  flow.  As  pointed  out  in 
section  1,  the  mixing  layer  exhibits  properties  that  are  independent  of 
downstream  location  when  scaled  with  the  appropriate  local  length  scale. 
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4.  Conclusions  and  Recommendations 

A  thorough  measurement  of  the  structure  and  dynamics  of  the  streamwise 
velocity  in  a  turbulent  plane  mixing  layer  was  taken,  using  a  single-sensor 
hot-wire  probe.  The  results  of  the  measurements  were  compared  against  a 
theoretical  computation  of  a  similar  flow  by  Rogers  and  Moser  (1994).  This 
investigation  demonstrated  that  the  wind  tunnel  facility,  with  the  splitter 
plate  insert  and  screen,  can  create  the  turbulent  plane  mixing  layer  envi¬ 
ronment  with  standard  velocity  statistical  properties  similar  to  that  in  the 
numerical  experiment  by  Rogers  and  Moser  (1994).  I  anticipate  that  sta¬ 
tions  further  downstream  of  the  last  measurement  station  will  be  well  into 
the  self-similar  region  described  by  Rogers  and  Moser. 

These  results  dictate  that  the  use  of  multisensor  probes  with  4, 9,  or  12  sen¬ 
sor  arrays  to  measure  the  velocity  vector  components,  the  velocity  gradi¬ 
ents,  and  the  vorticity  would  be  of  high  quality  and  value.  These  probes 
are  capable  of  measuring  minute  volumes,  such  that  the  probe  dimensions 
are  approximately  of  the  same  order  as  the  anticipated  Kolmogoroff  length 
scale  (using  the  isotropic  turbulence  assumption  for  calculation  of  the  dis¬ 
sipation  rate).  Thus,  high-quality  measurements  could  be  made  of  the  tur¬ 
bulence  properties  and  associated  statistics.  Conditional  averaging  could 
additionally  be  used  to  evaluate  the  distributions  of  turbulent  kinetic  en¬ 
ergy,  enstrophy,  dissipation  rate,  amongst  other  properties.  Such  measure¬ 
ments  could  be  compared  against  the  results  pf  prior  literature,  and  could 
then  be  extended.  Such  extensions  would  enhance  the  present  understand¬ 
ing  of  turbulence  in  coherent  structures,  increasing  the  capability  of  im¬ 
proving  turbulence  models  for  computational  fluid  dynamics,  and  would 
contribute  to  the  body  of  knowledge  in  turbulent  fluid  dynamics. 
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